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ABSTRACT 



Context. We use optical photometric and polarimetric data of BL Lacertae that cover a period of 22 years to study the variability of 
the source. 

Aims. The long-term observations are employed for establishing parameters of a stochastic model consisting of the radiation from a 
steady polarized source and a number of variable components with different polarization parameters, proposed by the authors early. 
Methods. We infer parameters of the model from the observations using numerical simulations based on a Monte Carlo method, with 
values of each model parameter selected from a Gaussian distribution. We determine the best set of model parameters by comparing 
model distributions to the observational ones using the x~ criterion. 

Results. We show that the observed photometric and polarimetric variability can be explained within a model with a steady source 
of high polarization, ~40%, and with direction of polarization parallel to the parsec scale jet, along with 10±5 sources of variable 
polarization. 
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1. Introduction 

BL Lacertae is an active galactic nucleus (z=0.069) that is 
the prototype of a class of AGN, BL Lac's objects, which 
are characterized by strong dominance of non-thermal con- 
tinuum, extreme optical variability, and high optical po- 
larization ( [Angel & Stockman 1980| l. BL Lac objects belong 
also to the class of blazars with highly variable emission 
from radio wavelengths to y-rays and superluminal appar- 
ent velocities in radio jets. Recent theoretical and observa- 



tional results (e.g., [Vlahakis & Konigl 2004| [McKinney 2006 
[Marscher et al., 2008| l suggest that the jets are magnetically 
dominated, with a helical structure of the magnetic field close 
to the central engine and accretion disk where the jet forma- 
tion and acceleration take place. However, at some distance 
from the central engine the jets become kinetically dominated, 
with a turbulent magnetic field that is subject to ordering by 
shocks. If the non-thermal optical emission in blazars origi- 
nates close to the millimeter-wave VLBI core, as suggested 
bv lLister & Smidi 2000 Ga buzda et al. 2006IIJorstad et al. 2007) 
then the optical polarization behavior should reflect the dynam- 
ics of the magnetic field in this region. 

On long-term timescales, many BL Lac objects show a pref- 



erential direction of optical polarization (Hagen-Thom 1980 



Hagen-Thorn et al. 2002 - Paper 1) closely aligned with the di- 
rection of the parsec scale jet jJorstad et al. 20071 . as well as 
some coiTelations between total intensity and polarization pa- 
rameters (Paper 1). The data on OJ 287 (Hagen-Thorn 1980) and 
BL Lac (Paper 1) permit to suggest a phenomenological model 
of variability for BL Lac object: there is a constantly acting 
source of polarized radiation responsible for the preferential di- 
rection of observed polarization and a number of sources with 
randomly distributed polarization parameters and intensities that 
are superimposed on the steady source. 
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In this paper we use simultaneous photometric and polari- 
metric measurements of BL Lac obtained over 22 years (for de- 
tails see Paper 1) and develop a method to derive the optimal 
number of the superimposed sources as well as the properties of 
the steady and superimposed sources that give the best fit to the 
data. We apply the results to describe the magnetic field structure 
in the optical emission region. 

We recall that the relative Stokes parameters and p,, are 
directly measured from polarimetric observations, and they are 
used to derive the degree of polarization P - (p^+Py)^^^ and po- 
sition angle of polarization ©o = 0.5 arctan(/?v/ p.^). Photometric 
measurements were performed independently and were cor- 
rected for contribution of the host galaxy. The transformation 
from magnitudes m to intensities / was performed with the ab- 
solute calibration of |Mead et al. (1990)| Polarimetric and pho- 
tometric measurements obtained on the same night were used 
to compute the absolute Stokes parameters I,Q - Pxl, and 

U = Pyl. 

2. Stochastic IVIodel 

Figures 1-6 {left column) show the observed dependences of po- 
larization parameters on Stokes parameter /. Figure 7a presents 
the plot between the degree of polarization and position angle of 
polarization, while Figure 8a gives the dependence between the 
relative Stokes parameters. 

The aim of this work is to reproduce all dependences ob- 
tained from the observations with a model consisting of a steady 
source of polarized radiation with parameters /'',/7j,pJ and n 
variable sources with parameters I],plj,plp where i - 
At a given date the observed optical emission is represented by 
Stokes parameters Imod, Qmod, and Umod that are computed as the 
superposition of the constantly acting and variable sources (see 
Appendix for the basic equations). The total number of dates is 
451. 
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All parameters are assumed to be normally distributed 
around their mean value (mathematical expectation M) with 
standard deviation cr. There are three constraints on the param- 
eters of the steady source: 1) its flux must be less than the min- 
imum observed flux, 2) its polarization direction must yield the 
preferred direction of polarization, which is Q'^^"^ - 24°, and 3) 
its degree of polarization must be sufficiently high to ensure a 
high observed polarization al low flux levels. 

Determination of optimal model parameters was accom- 
plished by a Monte Carlo method. After multiple realizations 
of the model with different sets of parameters, we accept a set of 
parameters that produces the best agreement with the observed 
dependences as the most probable model. The initial criterion 
for rejecting obviously inadequate models and outlining possible 
sets of parameters was the visual similarity between the modeled 
and observed correlations. 

The details of modelling were as follows: 1) the 
space of Stokes parameters {/, Q, U] was separated 
into K non-intersecting regions, each containing about 
10 points from the observed set of 451 points; 2) 
from possible sets of parameters a realization of values 
M{n), <T{n), M(n, criF), M{pl), <T{pl\ M(p;), cr(p';), F, p% p^ 
was selected by chance (the value of p'^^, was calculated from the 
condition 0q = 24°, and it was assumed that cr(/7'i) = cr(p'j)y, 
3) using this realization we found by random chance the values 
n,r ,p\,pl„r ,p'^^,p'[, for 451 cases; 4) using equations (lAll i - 
( |A7| l for each of the 451 cases, we determined Im„d, Qmod, Umod\ 
and 5) we compared the number of points in each of K regions 
derived in the model, MODj (i - 1, . . . , K), with the number of 
observed points, OBS ,■ (/ = 1 , . . . , K), by calculating the value 



X 



(MODj -OBS if 
OBSi 



(1) 



After repeating the third step 1000 times, we have found 
for the realization chosen in the second step. Regarding 
as a function of parameters of the second step, we have deter- 
mined the set of parameters that gives a global minimum of x^ 
by the genetic algorithm of optimization (Goldberg 19891. We 



used the PyGene library to implement the algorithm. Beginning 
with the second step, the calculations were performed 200 times. 
The whole algorithm was repeated for 3 different separations of 
{/, Q, f/)-space. This gives 600 values of parameters. The mean 
values that give the most probable realization are as follows: 
M(n) = 10 + l,o-(n) = 5 + 1,M(/') = 0.12 + 0.04,cr(/0 = 
0.15+0.05, M(p';) = 2+l,cr(p';) = 16+3, M(/:'P = 2±l,cr(p;;) = 
18 + 3,/" = 0.19 + 0.03,pJ. = 27 +5,p^ = 28+4. 



cr(/'') is higher than M{P') in some cases we formally may ob- 
tain by chance a negative value of the intensity /''; we exclude 
such cases. 

The proposed stochastic model agrees very well with the 
phenomenological model suggested in Paper 1 . The model gives 
a high degree of polarization for the steady source, P ~ 40%, 
with direction of polarization along the jet. This can be pro- 
duced by a toroidal component of a helical magnetic field in the 
magnetically dominated part of the jet. Such a scenario implies 
that the optical emission most likely comes from the accelera- 
tion zone of the jet, as suggested by |Marscher et al., 2008| who 
recently have observed in BL Lac steady rotation of the position 
angle of optical polarization by about 240° over a five-day inter- 
val. The rotation can be explained within a helical field structure 
located in the acceleration zone -0.4 pc upstream of the mm- 
wave VLBI core. However, the intensity of the steady source is 
only ~50% higher than the most probable intensity of a random 
variable source, while at a given moment we expect ~10 such 
sources to operate simultaneously. The most probable polariza- 
tion of variable sources ~3%, which implies that each variable 
source is an emission region with a significant number of cells 
with randomly oriented magnetic field. This means that, in gen- 
eral, the magnetic field is dominated by a turbulent component 
that masks the toroidal component. However, crp^j. and crpy are 
fairly high, which requires ordering of the magnetic field in the 
emission regions of the variable sources at some times. The lat- 
ter can be connected with propagation of disturbances along the 
jet and formation of shocks in the optical emission region. 



4. Appendix 

Because of additive nature of absolute Stokes parameters, the 
main equations for obtaining the model values are as follows: 



^mod 



= r + 



n 

!=1 

From ( lAll i - ( IA3l l for relative Stokes parameters we have 



Pxmod ~ 



(Al) 



(A2) 



(A3) 



(A4) 



3. Results and Conclusions 

We have constructed the dependences between the polarization 
parameters and intensity similar to the observed dependences 
using such set of parameters: M(n) - 10, cr(«) - 5,M(P') - 

0.12, cr(/') = Q.15, M{p\) = 2,cr(/?';) = i6,M(p;;) = 2,cr(;?;;) = 

18,/'' = 0.19, = 27,/?;: = 28. Figures lb-8b give the modeled 
dependences obtained for one realization of parameters taken at 
random from this set. There is an excellent agreement between 
the observed and modelled dependences. The similarity is ex- 
cellent for other realizations of this set of parameters as well. 
Figures I9I10I as examples, give two dependences obtained for 
two more different realizations. In addition, since the value of 



Py mod 



pr + i:up;,ij 



p + n,ij 

and obviously, 



Pmod 



}od'' 



0, 



Vmod 



■ arctan 



Py mod 
Pxmod 



(A5) 



(A6) 



(A7) 
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;s. 1-7. Observed dependences (left column) and modeling results (right column) between polarization parameters and intensity. 
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Fig. 8. Relative Stokes parameters: observed (left) and modeled (right). Mean values are indicated by open squares. 
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Fig. 9. The dependence of the degree of polarization on flux for 
two different realizations. 
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Fig. 10. The dependence of the degree of polarization on its di- 
rection for two different reaUzations. 



